Quantum Brachistochrone 
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We present a general framework for finding the time-optimal evolution and the optimal Hamil- 
tonian for a quantum system with a given set of initial and final states. Our formulation is based 
on the variational principle and is analogous to that for the brachistochrone in classical mechanics. 
We reduce the problem to a formal equation for the Hamiltonian which depends on certain con- 
straint functions specifying the range of available Hamiltonians. For some simple examples of the 
constraints, we explicitly find the optimal solutions. 
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In quantum mechanics one can change a given state 
to another by applying a suitable Hamiltonian on the 
system. In many situations, e.g. quantum computation, 
it is desirable to know the pathway in the shortest time. 

In this paper we consider the problem of finding the 
time-optimal path for the evolution of a pure quantum 
state and the optimal driving Hamiltonian. Recently, a 
growing number of works related to this topic have ap- 
peared. For instance, Alvarez and Gomez [lj showed that 
the quantum state in Grover's algorithm |2(, known as the 
optimal quantum search algorithm Q , actually follows a 
geodesic curve derived from the Fubini-Study metric in 
the projective space. Khaneja et al. and Zhang et 
al. y| , using a Cartan decomposition scheme for unitary 
operations, discussed the time optimal way to realize a 
two-qubit universal unitary gate under the condition that 
one-qubit operations can be performed in an arbitrar- 
ily short time. On the other hand, Tanimura et al. 
gave an adiabatic solution to the optimal control problem 
in holonomic quantum computation, in which a desired 
unitary gate is generated as the holonomy corresponding 
to the minimal length loop in the space of control pa- 
rameters for the Hamiltonian. Schulte-Herbruggen et al. 
Jl exploited the differential geometry of the projective 
unitary group to give the tightest known upper bounds 
on the actual time complexity of some basic modules of 
quantum algorithms. More recently, Nielsen Q intro- 
duced a lower bound on the size of the quantum circuit 
necessary to realize a given unitary operator based on 
the geodesic distance, with a suitable metric, between 
the unitary and the identity operators. However, a gen- 
eral method for generating the time optimal Hamiltonian 
evolution of quantum states was not known until now. 

Here we are going to study this problem by exploit- 
ing the analogy with the so-called brachistochrone prob- 
lem in classical mechanics and the elementary properties 
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of quantum mechanics. In ordinary quantum mechan- 
ics the initial state and the Hamiltonian of a physical 
system are given and one has to find the final state us- 
ing the Schrodinger equation. In our work we generalize 
this framework so as to optimize a certain cost functional 
with respect to the Hamiltonian as well as the quantum 
states. The cost functional quantifies the efficiency to get 
the target state from a given initial state and depends on 
the physical situation. In this paper we focus our atten- 
tion on the time optimality but it is straightforward to 
generalize our methods to other cost functions. 

In the brachistochrone problem one has to find the 
shape of a friction-free tube connecting two points and 
with a particle running inside subject only to homoge- 
neous gravity. The solution, a segment of a cycloid, can 
be found using the variational principle for the evolution 
time T[a;(£)] = J — , where the parameter s specifies the 
length of the tube from the initial to the current position 
x{t) of the particle, i.e. ds 2 — \dx\ 2 . The magnitude 
of the particle's velocity is v :— = \J%E — V(x))/m, 
where E is the conserved energy and V is the gravita- 
tional potential. 

Let us now move to the time optimization problem 
in the quantum case. We want to minimize the total 
amount of time necessary for changing a given initial 
state \tpi) (belonging to an n-dimensional Hilbert space 
Tt) to a given final state \ipf), by suitable choice of a (pos- 
sibly time-dependent) Hamiltonian H(t). In our problem 
the quantum state \tp(t)) and the Hamiltonian H(t) are 
the dynamical variables, and the action is defined as 



S(ip,H,(/>,\) = / dt 



(m-p)m 



AE 

i(<j>\ip) + ((t>\H\ip) + c.c. 
\(ttH 2 /2 - uj 2 



(1) 



Here P(t) := \%p(t))(%p(t)\ is the projection to the state 
\ip(t)}, H := H — (TrH)/n is the traceless part of the 
Hamiltonian, (AE) 2 := (i(j\H 2 \if)) - {i^\H\^) 2 is the en- 
ergy variance and uj is a given nonzero constant (an over- 
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dot denotes differentiation with respect to the time t). 
We have chosen units in which Planck's constant H is 
equal to one. 

The first term in the action gives the time dura- 
tion for the evolution of \ip(t)), expressed in terms of the 
Fubini-Study line element ds 2 = (di/j\(l — P)\dip) on the 
projective space CP™ -1 . The second term guarantees, 
through the Lagrange multiplier \(f>(t)) G P, that \tp(t)} 
and H{t) satisfy the Schrodinger equation and that the 
squared norm = 1 is conserved. The third term, 

through the Lagrange multiplier A, generates a constraint 
for the Hamiltonian. Such a constraint is necessary be- 
cause otherwise one would be able to find a path with ar- 
bitrarily small time duration just by rescaling the Hamil- 
tonian as H i ► aH, with a > 1, to make the energy 
fluctuations AP large. This corresponds to the fact that 
physically only a finite amount of resources (e.g., a finite 
magnetic field) is available. Here we consider a typical 
example which we call the isotropic constraint (we will 
consider more general constraints later). The constraint 
is imposed on P rather than P, since the difference be- 
tween the highest and the lowest energy levels in P, and 
not the value of the energy levels themselves, is important 
for the physical system. The problem should be math- 
ematically formulated on the projective space CP™ -1 
rather than on P, because the overall phase of the state 
|^) is of no significance in quantum mechanics. In fact, 
the action is invariant under the U(l) gauge trans- 
formation (\ijj),H,\(f)),\) i-> (e~ l0 \ij),H + 6,e~ ie \(l)),\), 
where 0(t) is a real function. Note that the Hamiltonian 
P plays the role of the gauge potential (there is also an- 
other symmetry, \<j)) t— > \<j)) + where 7 is a real 
constant). 

Let us now derive the equations of motion. The vari- 
ation of Q with respect to ((f>\ leads to the Schrodinger 
equation 



(2) 



In particular, this implies (^|(1 — P)|?A) = (AP) 2 , or 

ds = AE dt, (3) 

which was found by Aharonov and Anandan [tj and leads 
to a rigorous formulation of the time-energy uncertainty 
principle. The variation with respect to A gives the con- 
straint TrP 2 = 2ui 2 . The variation with respect to 
upon using J2J), yields 



\ 2(AP) 2 



i\6)+H\6)=0. 



(4) 



Finally, the variation with respect to P, after use of 
implies 



{P,P}-2(P)P 
2(AP) 2 



AP + MM) =0 (5) 



where a bracket (•) denotes the expectation value with 
respect to and {A, B} = AB + BA. Equation (gj and 



the trace of (0 imply that (ip\<j>) is a purely imaginary 
constant. Then, the expectation value of 10 gives (H) = 
0, which is equivalent to (P) = (Tr H)/n, or H = H — 
(H). Applying (J5J to |^), we have 



/ 1 



-A )H + {i>\4>) 



\2(AE) 2 

and inserting © back into ijHJ), we obtain 
H = HP + PH. 



(6) 



(7) 



Furthermore, the energy variance is constant, i.e. 



(AP) 2 



(P 2 ) = TrP 2 /2 



Substituting (JHJ into 



Q, we have (AP)' \ip) — which, after multiplication by 
(4>\H, implies that A is constant. We then obtain 



HI 



0, 



(8) 



where we have introduced := exp[i J* dt{H)]\ip). In 
conclusion, the equations to be solved have reduced to 
and ||8J|. Equation J7J gives an expression for the 
optimal Hamiltonian and JSJ gives the optimal time evo- 
lution of the quantum state. In fact, the state \ip) satisfies 
the Schrodinger equation with Hamiltonian H, and (JJJ) 
implies 



H = i{\^)^\-\i>){^\). 

The derivative \ip) is orthogonal to | 
Therefore, equation © reads 

(1-P)|J) =0, 



(9) 



because (P) = 0. 



(10) 



where P = IV'XV'I = P- This is the geodesic equation for 
the Fubini-Study metric on CP™ -1 , which is suggested 
by the observation that the first term in the action 
becomes / ds for constant AP. We also easily see that 

P = from © and lfTU|) . 



One can solve equation (|10|l using (AP)^ 



find- 



ing: 



\ip(t)) = cosut 1-0(0)) 



sin u)t 



(11) 



It is then easy to rewrite and H in terms of the 
Gram-Schmidt orthonormalized initial state \ipi) and fi- 
nal state \ijj'f) as 

\tj){t)) = cosLot\^i) + sin^|^), (12) 
P = w(|^)(^|-|^)(^|). (13) 

As a result, the whole Hamiltonian is given by H(t) = 
H + (H(t)), where (H(t)) is an arbitrary real function 
corresponding to the gauge degree of freedom. The opti- 
mal time is T = arccos \ (ipf\ipi)\. 
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Let us now generalize the quantum brachistochronc 
problem to the case of a more general set of m constraints 
for the Hamiltonian H(t). This extension is physically 
relevant since, for example, interactions between more 
than two qubits are not controllable in practical cases. 
Instead of Q we consider the action 



S{i/),H,\,n) = / dt 



y/(m-P)\ij>) 
AE 



m 

+ {4>\H\i>) + CC.) + E A< Va(#) 



(14) 



a magnetic field with x and y components only. By i|18|) 
and Tr F = 0, we have 



F = Ai# + \ 2 <r z = UF(0)UK 



(22) 



Eliminating H in 122|) via if/ = i?Z7 and solving (|22p. we 
get J7 = cxp[ixcr 2 ] exp[— iP(0) J d£/Ai], where x(t) := 

/oCWAi)di. From l(2T!l) - (|2*2^l and the above formula for 
£7, we find that the Xj are constants and we can simplify 
U and H as 



where the f a (a — 1, ...to) are functions mapping a Her- 
mitian operator into a real number. As we have already 
discussed, the most natural case to consider is when the 
f a are actually functions of H, but the argument below is 
valid also for arbitrary f a . Among the equations of mo- 
tion, the constraint now generalizes to f a (H) — 0, while 
and JHJ) become 



F = FP + PF, 
F + i[H, F]) =0, 



with the operator 
F(H): 



\5H "5H' 



P 



(15) 
(16) 



(17) 



In particular, equation (|15fl guarantees that TrP = 
(F) = 0. Note that, for m = 1 and / = TrP 2 /2 - u? , 
equations H15|) and l|16[) reproduce (JJJ and JSJ) of the 
isotropic case. 

We can formally integrate (|15|) and 116(1 to obtain 



F = UF(G)U\ 



(18) 



where F(0) is a constant Hcrmitian operator which sat- 
isfies F(0) = {F(0), P(0)} and the unitary operator U is 
given as a functional of H , 



U[H](t) :=fe- l ti" d \ 



(19) 



with T the time-ordered product. Thus, given the con- 
straints fa(H) — 0, one can explicitly write the left hand 
side of (|18fl as a function of H and solve (|18(l to obtain 
the optimal Hamiltonian H . This is our main result. 

As an explicit example of this general framework we 
may consider the following case of one qubit subject to 
the two constraints 



h{H) : = TyH 2 /2-u 2 
j 2 {H) : = Tr(Ha z ) = 0, 



0. 



(20) 
(21) 



where a.j (j — x, y, z) are the Pauli matrices. This cor- 
responds to the physical situation where one can apply 



jj _ e intaz e -i[H(o)+naz]t 



(23) 
(24) 



where :— A2/A1. 

Let us consider the case in which the initial state is on 
the equator of the Bloch sphere CP 1 , (a z ) = 0. With- 
out loss of generality, we may choose P(0) = (1 + tT x )/2, 
i.e. the spin initially points towards the positive x direc- 
tion. Then the condition F(0) = {F(0),P(0)} and the 
constraints (|20|l and l|21|) imply that P(0) = — uj(T y , and 
from 112:') ll and l|24|l we finally obtain 



(*->(*) 



cos 2Q,t cos 2Q't - 
- sin 2ttt cos 2fl't 



n T sin 2VLt s\n2Q't 



it 



^rCOs2ntsm2n't 
7 sin 20,'t 



H{t) = -<r-B(t) ; B(t) 




(25) 
(26) 



where ft' := ^/uj 2 + fl 2 . We can interpret B(i) in 
()26|l as a magnetic field rotating with angular veloc- 
ity 2Q,. Equation (|25|l is not a geodesic on CP 1 un- 
less f2 = 0, in which case the orbit is a great circle 
in the xz-plane. The energy fluctuation (which is also 
the speed of the state) is now time-dependent, AE(t) = 



l-(§sm2n't) 



1/2 



We notice that the constraint 

(|21ll in general reduces the speed of the state. 

If the final state \ipf) is given, the angular velocity 2S1 
and the time duration T are also determined. For in- 
stance, let us assume P(T) — (1 — <J X )/2 1 which is the 
antipodal point of P(0) in CP 1 . Then we have the con- 
ditions 2|fi|T = kn and 2il'T = In, where k and I are 
integers such that I > k > and k + I is odd. Thus we 
get 



\T=^l 2 -k 2 , 



VP 



(27) 
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FIG. 1: Locally optimal curves on the Bloch sphere. The 
solid, thick dotted and thin dotted curves correspond, respec- 
tively, to the solutions with \u>\T — ^, -| y/3 and f^5. They 
have zero, one and two nodes along the great circle on the 
a;y-plane, respectively. 

The concrete solutions for J57J), in ascending order of 

T, are (\u>\T, |S|) = (f , 0) , (f y/3, ^) , (f A ^) , • • • . 

These solutions respectively have zero, one, two • • • nodes 
along the great circle (a z ) — 0. The first one is a geodesic, 
while the others are not. On the other hand, if the final 
state is not in the xz-plane, there is no geodesic solution. 
At this stage we have to caution the reader that solutions 
may be only locally optimal, and we have to find out the 
globally optimal one by comparing time durations T . In 
the above case, the first curve without nodes is globally 
optimal. The second curve becomes globally optimal if 



the final state is its first node. The solutions are depicted 
in Fig. 1. 

In summary, in analogy to the classical brachistochronc 
problem we have formulated a variational principle to 
find the optimal Hamiltonian and the optimal quantum 
state evolution, for given initial and final states and a set 
of available Hamiltonians. As a particular application of 
our methods, one might want to first evaluate the optimal 
Hamiltonian, e.g. by means of a classical computer, and 
then perform the fastest possible quantum experiment or 
quantum computation. As a future development of the 
present research we do not see any obstacle to generalize 
our formulation to mixed states. The relation to gate 
complexity in the standard paradigm of quantum com- 
putation remains to be investigated, though our point 
of view is that the time complexity (see, e.g., Schulte- 
Herbriiggen et al. J\ and references therein) is more 
physical or even more practical. For example, equations 
(I12H and (|13|) essentially solve Grover' s search problem 
QlQl- However, to show the square-root speed up with 
respect to the classical case one still has to identify the 
computational step corresponding to a single oracle call. 
We hope that the present variational approach to the 
time optimality problem in quantum mechanics opens 
up novel systematic investigations of optimal quantum 
computation. 
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